Semi-Meissner state and neither type-I nor type-II superconductivity in 

multicomponent systems. 



Egor Babaev^'^ and Martin Speight^ 
^Laboratory of Atomic and Solid State Physics, Cornell University Ithaca, NY 14853-2501 USA 
'^Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 
^ School of Mathematics, University of Leeds, Leeds LS2 9JT, UK 

Traditionally, superconductors are categorized as type-I or type-II. Type-I superconductors sup- 
port only Meissner and normal states, while type-II superconductors form magnetic vortices in 
sufiiciently strong applied magnetic fields. Recently there has been much interest in superconduct- 
ing systems with several species of condensates, in fields ranging from Condensed Matter to High 
Energy Physics. Here we show that the type-I/type-II classification is insufficient for such multi- 
component superconductors. We obtain solutions representing thermodynamically stable vortices 
with properties falling outside the usual type-I/type-II dichotomy, in that they have the following 
features: (i) Pippard electrodynamics, (ii) interaction potential with long-range attractive and short- 
range repulsive parts, (iii) for an n-quantum vortex, a non-monotonic ratio E{n)/n where E{n) is 
the energy per unit length, (iv) energetic preference for non-axisymmetric vortex states, "vortex 
molecules". Consequently, these superconductors exhibit an emerging first order transition into a 
"semi-Meissner" state, an inhomogeneous state comprising a mixture of domains of two-component 
Meissner state and vortex clusters. 



The formation of vortices in type-II superconductors 
subjected to a magnetic field is one of the most remark- 
able phenomena occuring in condensed matter. In all 
type-II superconductors (i.e. superconductors where the 
GL parameter, which is the ratio of the magentic field 

Binetration length to the coherence lenth, is k > l/\/2 
) these vortices share a set of properties: an n-quantum 
vortex is unstable with respect to decay into n one- 
quantum vortices, two vortices have purely repulsive in- 
teraction, and invasion of vortices under normal condi- 
tions is a second order phase transition characterised by 
a critical value of the external magnetic field Hd- 

Vortices as solutions of the GL equations also exist for- 
mally in a type-I superconductor {k < l/\/2), but these 
vortices are thermodynamically unstable. The special 
case K = 1/ V2 is also very interesting since, at this value 
of K, vortices do not interact [3, Q. One should note, 
however, that in real life systems the situation is more 
complicated; experiments % show that in certain mate- 
rials with k « 1/^/2 there might exist a tiny attractive 
force between vortices at a certain distance. Such an in- 
teraction was reproduced in a modified one-component 
GL model with additional terms in the regime k « 1/ \/2 
[3- We should also mention a long range van der Waals 
- type vortex attraction in layered systems produced by 
thermal fiuctuations or disorder 0. 

Besides superconductivity, the vortex concept has a 
direct counterpart in High Energy Physics, called the 
Nielsen-Olesen string. Such strings have been considered 
in cosmology where they are expected to form during 
a symmetry breaking phase transition in the early uni- 
verse. There also exists a similar type-I/type-II division 
of semilocal cosmic strings in the Higgs doublet model 

i 

Recently, multicomponent superconducting systems 



have attracted increasing interest in areas ranging 
through metallic superconductors, hydrogen in extreme 
conditions and color superconductivity in dense QCD 
0, 0, . Below we consider a generic two-component 
superconductor (TCS), showing that it allows a novel 
type of thermodynamically stable vortices, whose elec- 
trodynamics is of Pippard type with respect to one of 
the order parameters and which have non-monotonic in- 
teraction energy for a wide range of parameters as an 
intrinsic feature. We will show that, as a result of this, 
such a TCS displays very unconventional magnetic prop- 
erties which have no counterparts in single-component 
systems, and do not fall into the standard type-I/type-II 
classification. 

In the simplest case, the TCS (related to the two- Higgs 
model (12i ) can be described by the following GL energy 
density: 
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where = |*Q|e"P° and 

Ea=i,2-^a|*"l' + f In the' model Q, 

vortices with phase winding in only one field have 
logarithmically divergent energy per unit length if ?/ = 
and linearly divergent if 77 ^ |^. Here we do not 
consider the effect of thermal fluctuations, so we can, 
without loss of generality, restrict attention to vortices 
with ipi — 4>2 = const in the case 77 = 0. These have 
finite-energy per unit length. Moving 77 from zero will 
merely change the core sizes of such vortices, so this is 
a trivial extension. The same applies to other possible 
additional potential terms in l(T|l. 

The model Q possesses three characteristic length 



2 



scales: two coherence lengths £,a = h/y/Aniaha 
and the magnetic field penetration length A = 
(c /V8^e )[|^i|Vmi + 1*217^2]"'/', where = 
^/ha/ca- It will also be convenient to define length 
parameters Aq. = c^JnTZ/ {-s/STreY^ a\), and thermody- 
namic critical magnetic fields for the individual con- 
densates, Hct(a) = ^o/i'^V^'n'S.a^a)- When the indi- 
vidual condensate is of type-II, we denote its first 
and the second critical magnetic fields as -ffci(i) = 
$o/(47rA?)[log(Ai/a) +0.08] and ff,2(i) = $o/(2^e?), 
where $o is the magnetic flux quantum. 

Let us consider the magnetic properties of a TCS in 
several regimes. In the simplest cases, when ~ ^2 >> 
A and when $1 ~ ^2 << A, the magnetic properties of 
the TCS parallel those of single-component type-I and 
type-II superconductors correspondingly. However, we 
find that if one of the condensates is of type-I while the 
other is of type-II, the TCS in an external field has a 
much richer phase diagram. 

In the case when ^1 << Ai and ^2 >> A2, and -ffct(2)j 
the thermodynamic critical magnetic field for the type-I 
condensate ^'2, is much higher than H^2(i) for the type-II 
condensate the system undergoes a first order tran- 
sition from a normal state immediately into a TCS state. 
This is because, at fields higher than 7?j,t(2)j the system 
does not allow nontrivial solutions of the linearized GL 
equation for ^'i, while, when the field is lowered below 
Hct{2) J there appears a transition immediately into a TCS 
state. In this state the magnetic field is screened in the 
bulk of the sample largely due to surface current ^'2, and 
nothing can preclude the appearance of the second con- 
densate ^'i in the bulk of the sample, even if the applied 
field is much larger than -ffc2(i)- 

Now consider the regime occuring when Ai/^i > 
l/V^;X2/^2 < l/\/2;^2 > A. Then a vortex solution 
should have an extended core associated with the con- 
densate ^'2 which exeeds the penetration length (the 
crucial question of whether such vortices can be ther- 
modynamically stable will be answered below). For a 
TCS, the vortex energy consists of the energies of the 
cores, the kinetic energy of the Meissner current and the 
magnetic field energy. The magnetic field energy and 
the kinetic energy of the screening current are given by: 
F„, = (I/Stt) J d^xU^ + (I/Stt) J d^xXlgicurm.)^. Here 
we stress that, in the present case, |4'2(x)p varies slowly 
over the London penetration length A. The magnetic 
field is screened at a distance from the core which is 
smaller than ^2. This means that only a depleted den- 
sity of Cooper pairs of the condensate ^'2 participates 
in the screening of the magnetic field. Thus, one can- 
not use the London penetration length A, but should in- 
troduce an effective penetration length AeS = [1/A^ + 
l/A2(x)]-i/2 where_A2(x) = c^/(V8^e|*2(x)|) > 
A2 = Cy/rfi2/ (\/87re|^'2|). In the case when Ai << ■^2 and 
Ai is much smaller than the Pippard length of ^'2 we have 



Aoff ~ Ai which corresponds to the situation when the 
magnetic field is screened mostly by condensate On 
the other hand, in the case when ^2 >> Ai >> A2 >> ^1 
the magnetic field can be screened at the scale of the 
Pippard penetration length A = Af « {'^2^2)^^^ of the 
condensate 5*2 (this expression is valid when Af << Ai). 
In the above expression for Fm, we cut off integrals at 
the distance ^1 from the center of the core in order to 
obtain an estimate of the vortex energy with logarithmic 
accuracy. Then the energy per unit length of a one-fiux- 
quantum vortex is 

where Vca are the energies of the cores per unit length 
which are of order of magnitude of [core size] x 
[condensation energy] . The estimate of the core energy 
can also be expressed as: 

8^ " 8^ c2m„ ~ A{A7rJ XI ^ > 

Consequently, the energy of two cores is Vd + Vc2 ~ 
(<I'o/87rA)2. Let us now assume that such vortices are 
thermodynamically stable (below we demonstrate nu- 
merically that this is indeed the case). Then a straight- 
forward calculation of the field H^i at which it becomes 
energetically favorable to let a single vortex into the su- 
perconductor gives: 
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However this characteristic field strength cannot be inter- 
preted as the first critical magnetic field as in type-II su- 
perconductors. In fact, in contrast to the type-II regime, 
invasion of vortices into a superconductor in the regime 
in question should be accompanied by a magnetization 
jump and be of first order. This is because the vortex 
has an extended core at length scale ^2 > AeS which gives 
rise to attraction between vortices U oc Kf){^/2r / £^2) for 
r > Aoff, where Kq, is the Bessel function. The attrac- 
tion originates in winning in condensation energy in the 
condensate ^'2 when outer cores overlap. As we shall 
see below, the interaction potential has a repulsive part 
at a shorter length scale of the order of Acff. Conse- 
quently, a lattice of vortices with a spacing determined 
by the minimum of the interaction potential is preferred 
over a system of widely separated vortices. So the en- 
ergy of a system of n vortices in this regime is mini- 
mized when vortices spontaneously form lattice clusters 
with overlapping outer cores. We should observe that 
H^i in the present situation is larger than the thermo- 
dynamic critical magnetic field of the type-I condensate, 
Hct{2) = <I'o/(2-\/27r^2A2). We stress that this does not 
mean that at the field l@J the condensate ^2 is com- 
pletely depleted, however. This is because, when the 
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applied field is close to H^^, the field is mostly screened 
by the supercurrent of the condensate which circu- 
lates along the sample's edge. Thus the vortex system is 
dilute, or, more precisely, the intervortex distance is only 
determined by the effective attraction. So in contrast 
to the usual type-I / type-II behavior, the TCS in this 
regime displays a first order transition into an inhomo- 
geneous state consisting of clusters of vortices, where the 
order parameter '^2 is depleted due to the overlap of outer 
cores. So superconductivity in these vortex "droplets" is 
dominated by the order parameter and is essentially a 
one-component superconductivity. Since the vortex den- 
sity depends on the applied field, while the intervortex 
distance is determined by the nonmonotonic interaction 
potential, there should be present, besides these clusters 
of vortices, domains of two-component superconductivity 
in the vortexless Meissner state. We call this the semi- 
Meissner state. The transition into the semi-Meissner 
state may be viewed as a vortex matter analog of the 
condensation of water droplets or a sublimation process 
in classical physics, with the external field playing the 
role of "pressure" . At a higher value of external field the 
system will transition from the semi-Meissner state to 
the regular Abrikosov lattice. This transition might be 
rather complicated because of the features of the vortex 
interaction potential. That is, the standard argument 
in favour of triangular lattice symmetry no longer holds. 
Different lattice symmetries for a given density of vor- 
tices are characterised by different numbers of neighbors 
and different nearest-neighbour distances, so one might 
construct a sequence of transitions between lattices of dif- 
ferent symmetries as a function of external field strength. 
Again an analogy between the external field strength and 
the role of pressure in classical physics might be invoked. 

The key question regarding the existence of this tran- 
sition and the semi-Meissner state is whether vortices in 
this regime are stable not only topologically but also ther- 
modynamicaly. Below we find an answer to this numeri- 
cally. 

For the purposes of numerical simulation of this 
system, it is convenient to rescale the fields. Let 
A = (2e/ftc)A, = e/cA/(87r/m„)5'a. Fur- 

ther, we introduce some new parameters = 

(27re^ c^/ft^)m^CQ, = (Aba/T^e^c^maCa), so that 

2 

the effective potential becomes Va = ^(m^ ~ IV'qP)^- 
Then the free energy density in the case 77 = is 

i^^^^ l|VxA|2+ ^ i|(v + zA)^j2^K.. (5) 

a=l,2 

In terms of the new parameters, the 3 natural length 
scales are the inverse masses of the photon and the 
two Higgs bosons, that is, \ = {u\ + u\)~^ , = 
(/iaUa)"^. This system supports radially symmetric so- 
lutions 

A = r-^a(r)(- sin 6*, cos 61), = cra(r)e'"^ (6) 



where a and Ua are real and satisfy the boundary condi- 
tions a(0) = (70,(0) = 0, a(oo) = — n, cTq, (00) = Wq. This 
ansatz reduces the field equations to an ODE system, 

a" - - {n + a){al + al) = 

< + ^-in + ar^ + ^aUul-al) = 0, (7) 

solutions of which may be found by means of a shoot- 
ing method similar to that used in (see remark 
Figures lA-lC were generated using this scheme, em- 
ploying a 4th order Runge-Kutta method with variable 
r step for the numerical integration, and with ro = 0.01, 
ri = 2, Too — 8. The regime of most interest is the one 
where ^i/V^ < A < C2/V2. To explore the dynamical 
issues of interest, we would like to make the disparity 
between ^1 and ^2 as extreme as possible. There is a 
numerical limit to the disparity of length scales we can 
achieve (see remark |l5j). but the phenomena of interest 
can be demonstrated within that limit. 

Figure lA shows the profile functions of a single 
two-component vortex in the neither type-I nor type-II 
regime, where the electrodynamics is of Pippard type in 
respect of component \['2. Note the disparity in healing 
lengths of the two condensate fields. Figure IB shows 
the energy of a single two-component vortex as a func- 
tion of ^2, with A and ^1, fixed, normalized by Eq, the 
energy of a one- component vortex with the same A and 
^ = ^1. Figure IC shows the energy per vortex of n = 1, 
n = 2 and rt = 3 co-centered vortices at fixed ^1, ^2, A, 
as a function of |\E'i|/|^'2|, normalized by Eq, the energy 
of a single one-component vortex with the same A and 
^ — ^1. When l^fil = |\1/2|, the n = 1 vortex is favored, 
but as |5'i|/|5'2| increases, first the n = 2 vortex then the 
n = 3 vortex become energetically favored. This is a very 
interesting property since such an effect has no counter- 
part in type-II or type-I one-gap superconductors (N=l 
Abelian Higgs model). We should note that the n = 2,3 
co-centered solutions are unstable with respect to forma- 
tion of non-axisymmetric n-quantum vortex molecules, 
as we shall now see. 

To compute the interaction energy of a vortex pair 
(shown in Fig. ID), one must go beyond the radially 
symmetric ansatz and resort to a lattice minimization 
method (see remark jl^)- Note that the intervortex force 
is attractive at long range but has a repulsive core, as pre- 
dicted. The only stable static two-vortex we find in this 
regime has nonzero vortex separation and broken axial 
symmetry. For an isolated one-quantum vortex, we find 
numerically that when ^1 <C A <C ^2 and Ai is much 
smaller than the Pippard length of the condensate ^2, 
the vortex energy per unit length tends asymptotically 
to the energy of a vortex in a single condensate '^i plus 
the core energy of the vortex of condensate ^'2. Noting 
that the thermodynamical critical magnetic field of ^'i 
is proportional to (CiAi)"^, while the core energy of the 
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FIG. 1: A: The numerical solution for the neither type-I nor 
type-II vortex. B: The energy of a double-core vortex as a 
function of the disparity of the two coherence lengths. C: 
The non-monotonic ratio of energy per unit length to number 
of flux quanta for axisymmetric n-quantum vortices. D: The 
non-monotonic interaction potential for double-core vortices. 



vortex in ^'2 is proportional to , this proves the ther- 
modynamical stability of the neither type-I nor type-II 
vortices whose electrodynamics is of Pippard type, and 
the existence of the semi-Meissner state. 

In conclusion, superconductivity in multicomponent 
systems has recently attracted much interest in the 
physics of condensed matter and beyond. Here we show 
that, in contrast to ordinary superconductors, multicom- 
ponent systems allow for thermodynamically stable vor- 
tices even in the Pippard regime. Moreover, their mag- 
netic properties in a certain range of parameters do not 
allow one to classify such a superconductor as type-II 
or type-I. Rather, it should legitimately be placed in a 
separate class. Such a novel type of superconductivity 
should be relevant for a variety of systems. For example, 
it is well known that disparity of coherence lengths oc- 
curs naturally in two-band superconductors. E.g. MgB2 
and Mgi_xAlxB2, while type-II superconductors not be- 
longing to the regime we consider in this paper, have sig- 
nificant disparity in coherence lengths Analogous 
situations might appear in mixtures of condensates with 
different pairing symmetries. Another candidate for this 
type of superconductivity is the projected liquid metal- 
lic state of hydrogen [iflj where this regime is expected 
to be realized under certain conditions. A similar situa- 
tion might also occur in the color superconducting state 
in quark matter ^ . Certain features of the considered 
state should be preserved and might be experimentally 
accessed for vortex stacks in layered systems when 
one layer is type-II and another is strongly type-I. 
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Since the system is singular at r = and r = 00, we 
shoot forwards from r = ro, small, and backwards from 
r — Too, large, applying a matching condition at an inter- 
mediate radius ri . The shooting parameters are the Tay- 
lor coefficients Si for the small r fields, and qi, the mag- 
netic dipole/scalar monopole charges of the large r fields 
fl3l : small r: a{r) = sor^ + • • • , cra(r) = Sar" + ■ ■ ■ 
large r: a{r) = -n + qorKi{^J u\ + u^t) + ■ ■ ■ ,Oa(r) = 
Ua — qaKoi/J-aUar) -|- ■ • • . The problem is thus reduced 
to zeroing an R'' valued function (the mismatch at ri 
of a, ai , (72, a', (jj , (J2) on R^, which may be solved by a 
Newton-Raphson method. 

When ^2 is large, the field ip2 has a long range, so ap- 
proaches its vacuum value only slowly. We must thus 
choose roo to be large. But then if is small, iti — ai (roo) 
is tiny, and large changes in qi axe required to control the 
shots. Hence, the Newton-Raphson scheme is close to sin- 
gular, and the numerical method becomes unstable. 
The idea is to fix a value of the vortex separation p, 
then minimize E = F subject to the constraint that 
ipc{-p/2,Q) = ipa.(p/2,Q) = 0, a = 1,2. We chose to 
place the system on a regular spatial lattice (spacing 
Sx = 5y = h) with forward differences replacing the 
partial derivatives occuring in F, then to solve the con- 
strained gradient fiow equation for E using the Euler 
method with time step 5t, until the configuration con- 
verged to a minimum (up to numerical tolerance, that is, 
until the energy loss per time step is less than Etoi). Fig- 
ure ID was produced using this scheme with /i = 0.1 and 
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grid to ensure that the ip2 field experiences no significant 
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